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Regular two-dimensional lattices of evanescently coupled waveguides may provide 
Ph ■ in the near future photonic components capable of combining interferometrically and 

, simultaneously a large number of telescopes, thus easing the imaging capabilities of 

optical interferometers. In this paper, the theoretical modeling of the so-called Dis- 
crete Beam Combiners (DBC) is described and compared to the conventional model 
used for photonic beam combiners for astronomical interferometry. The performance 
of DBCs as compared to an ideal ABCD beam combiner is discussed and applications 
to astronomical instrumentation analyzed. 

>: 
ov 

^t- 

00 
q 

06 
o 



2 



I. INTRODUCTION 



Long baseline, optical interferometry is the only sustainable option to achieve ultra high 
resolution imaging in astronomy [l?] ]. The retrieval of accurate interferometric images re- 
quires however a dense sampling of the spatial coherence function (complex visibility) of 
the light emitted by an astronomical target ^J, a task which can be eased by the simulta- 
neous combination of light collected by several telescopes. To this end, integrated optics 
can provide miniaturized multiple-telescope combiners which can be included in very com- 
pact assemblies, while offering highly stable an repeatable visibility measurements {3I, Isj]. 
State-of-the-art photonic beam combiners can deliver simultaneous combination of up to 
four telescopes and provide visibility measurements over 6 interferometric baselines The 
modular and compact design of such components has significantly reduced the complexity of 
instrumental design and relative calibration procedures, with considerable reduction of the 
time required to commissioning an instrument |Sj. Existing photonic beam combiners are 
fabricated with silica-on-silicon technology, which constrains the waveguides and couplers to 
a plane. This represents a significant complication of the design of components allowing the 
combination of a significantly larger number of telescopes, due to the necessary presence of 
cross-overs and increased sensitivity to fabrication defects(2|. In this respect, three dimen- 
sional (3D) photonic components, such as those fabricated by direct laser- writing technique 
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may provide a significant simplification o 



the design of interferometric beam combiners 
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22]. The degree of simplification in design 



through off-plane fiber connections/couplers 
and fabrication of 'on chip' astronomical interferometers can indeed be dramatic, as proved 
by the recent proposal by Minardi fc Pertsch Q of us,ng regular 2D-arrays of evanescent.y 
coupled waveguides (i.e. photonic lattices^) rather than a traditional cascade of couplers 
arranged in a 3D environment. The so called Discrete Beam Combiners (DBC), besides hav- 
ing a very simple design, have the potential of improving the scalability of integrated beam 
combiners to telescopic arrays of arbitrary size, while providing a slightly higher sensitivity 



respect to the ABCD 



20| or pupil remapping schemes jlOy. Moreover, the possibility to 



fabricate DBCs with direct laser writing allows the possibility of cost effective prototyping 
on a range of materials suitable for virtually any optical wavelength from UV to MIR. 

Aim of this paper is to provide an extended theoretical background of DBCs and a first 
analysis of the performance of the DBC scheme as compared to idealized beam combiners. 
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In the first section the conventional theoretical modeling of photonic combiners (Visibility 
to PhotoMetry, V2PM, [111]) is recalled and linked to a formalism specifically developed for 
the analysis of the DBCs (a-matrix approach). In Section 3, after a short recall of the DBC 
concept and the theory of evanescent coupling of optical waveguides, the method to derive 
the a-matrix for square arrays of waveguides is derived. In Section 4, the performance of 
ABCD combiners is compared to the DBC combining 3 to 6 telescopes. Section 5 concludes 
the paper indicating application perspectives of DBCs in astronomical interferometry. 



Interferometric beam combiners are used to encode phase variations into intensity vari- 
ations, suitable for the retrieval of the coherence properties of the combined optical fields. 
Ideally the output port of a beam combiner provides the interference signal of just one pair 
of the input fields combined by the device N input complex fields {A n }: 



Here, (•) indicates the time averaging, 3? is the real part, <pi 2 is the average phase difference 
between the fields and 7 and 7 are the amplitude and phase of the complex visibility of 
the interfering fields. More generally, the M measurables {I m } of an interferometric beam 
combiner are the real part of time averages of all possible products of N input complex fields 
{A n } with their conjugates (mutual coherences). This can be written in a compact form 



II. THEORY OF INTERFEROMETRIC BEAM COMBINATION 



A. The V2PM approach 



/ = (A X A\) + (A 2 A* 2 ) + 29* [(A.Al)} 
= (AiA\) + (A 2 A* 2 ) + 2 I Ax A\ | 7 cos(0i 2 + 7 ) 



(1) 
(2) 



using the V2PM (visibility to pixel matrix) formalism 11, 




N(N+1) 



Im = M{ E {V2PM} mn V n }. 



(3) 



n=l 
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Here, V2PM is a MxL complex-valued matrix, and the components of the vector V are all 
possible time averaged products of the input fields with their conjugates: 



V 



(A 1 A* 1 ) \ 
(AnA* n ) 

(AA* N ) 
(A 2 A*) 

{ (A N ^A* N ) ) 



(4) 



The elements of the V2PM matrix are calculated from the complex transmission coefficients 
of the input fields to the M outputs of the beam combiner. In existing planar photonic 
beam combiners exploiting the ABCD method Q], the matrix is formed by 4xL blocks, 
as apparent from Fig. 1 llj. The coherences are then usually obtained in the least square 
sense by pseudo-inversion methods. In a real beam combiner, the fabrication defects give 
raise to deviations from the ideal matrix of Fig. 1 which may take the form of i) amplitude 
variations of the non-zero elements, ii) deviation from the ideal phase of the complex matrix 
elements, and Hi) finite amplitude and phase of the zero elements (cross-talk). 



B. The quadrature approach 



It is possible to reformulate the mathematical model of the beam combiner in terms 
of products of real valued matrices and vectors. The method is based on considering the 
transformation of the real and imaginary parts (quadratures) of the complex field products 
which uniquely define the mutual coherence properties of the interfering fields. We can thus 
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FIG. 1. The V2PM matrix for an ABCD 3-beam combiner 



build a vector J of dimension N 2 which replaces the complex vector V: 



J = 



(AnA* n ) 

^(A!A* N ) 
M(A 2 A* 3 ) 

Q(A N ^A* N ) 
ZiA.A*) 

%(AiA* N ) 
Z(A 2 A* 3 ) 



(5) 
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FIG. 2. The ABCD 3-beam combiner in the a-matrix formalism 

Here, 9ft is the real part, and S the imaginary part of the complex field product. The output 
of the beam combiner is now described by the real matrix a of dimensions MxN 2 : 

N 2 

Ira ^ t ^miii (6) 
n=l 

For an ideal beam combiner giving 4 quadratures for the pair m n of the input fields, the 
corresponding block of matrix alpha is depicted in Figure 2. The change of the base of the 
coherences adds non-zero elements to each block. Again, the matrix has a non vanishing 
null space, thus the coherences can be obtained only by means of a pseudo inversion of the 
a matrix from the output intensity measurements. 



III. DISCRETE BEAM COMBINERS 



A. The concept 



As discussed in [12], two-dimensional regular arrays of NxN coupled waveguides can be 
used to determine uniquely the coherence properties of N fields {A n }. A conceptual scheme 
of the beam combination with arrays of waveguides is outlined in Figure 3 for the case N=4. 
In arrays of coupled waveguides, light injected into a waveguide spreads to neighboring sites 
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FIG. 3. A schematic view of the DBC suitable for interferometric combination of 4 telescopes. An 
array of 4x4=16 evanescently coupled waveguides is excited in the 4 points by the fields collected by 
the 4 telescopes. After suitable propagation in the array, a discrete excitation pattern of the array 
can be recorded (e.g. by 4x4 single pixel detectors glued to the end of the array of waveguides) 
which can be linearly related to the mutual coherence properties of the input fields through the 
inverse of the a- matrix (see text for details). 

upon propagation. By exciting several waveguides simultaneously, the field intensity of the 
excited waveguide mode as measured at the end of the array results from the interferometric 
mix of all input fields contributing in variable proportion, depending on the excited input 
site and length of the array. For a given configuration of input sites and array length, the 
discrete interference pattern can be related conveniently to the mutual coherences of all 
possible combinations of the input fields. 

As long as low index contrast fibers are concerned, modeling of light propagation can be 
carried out in terms of the coupled mode equations. By indicating with E n the complex 
amplitude at the peak of the mode propagating in n th waveguide, it is possible to describe 
the propagation of the fields along the longitudinal coordinate z by means of a system of 
coupled differential equations: 



dE, 



n 




(7) 



dz 



nm 



where the coupling coefficients c nm are proportional to the overlap integral of the normalized 
complex field transverse mode profiles u(x, y) of waveguides n and m. For an array of iden- 
tical waveguides with core index n core , substrate index n and effective propagation constant 



8 



/3, we can write: 

k 2 

_ I 2 _ 2\ ^0 
Cnm — l^corc n ) c^Q 

where the integral is evaluated over the cross section E of the core of the fiber. Strength of 
the inter-waveguide coupling can be varied by means of a suitable geometric arrangement 
of the fibers. Solutions of equation (7) can be obtained by direct integration or by means of 
the supermode decomposition of the system of equations. 



u D (x,y)u* m (x,y)dxdy, 



(8) 



B. Arrays of waveguides as combiners 

The a-matrix formalism is best suited to model the output of each waveguide for a 
given array of coupled waveguides, because it is possible to derive an exact expression of its 
elements from the solution of the coupled mode equation describing propagation of light in 
arrays of waveguides (Eq. 7). In fact, the peak intensity at the output of the m th waveguide 
of length z=L is given by: 



Im = (\E m \ 2 ) = ( 



M z 

fc=l 
N M 

= EE a m,/(j)< /(k ) (AjAl), (9) 

j=l k=l 

where and /(k) is a function mapping k=l...N onto the sites of the NxN waveguides array 
where the fields are coupled. The coefficients a n ,/(k) are the mode amplitudes at the 
end of waveguide n when a field of unit power is injected in site f(k) at z=0. They are 
function of the array geometry and sample length and are calculated by solving equation 
(1). The components of the a-matrix are then obtained straightforwardly from the following 
expressions: 

«n,k = K,/(fc)| 2 k = 1...3 (10) 
«n, P (j,k,M) = = 29?{a nJ(j) a* /(fc) } (11) 
tt n , q (j,k,M) = = -23{a ni/(j) a* )/(fe) }, (12) 

where the indices p and q are defined as: 

p(j, k, M) = j + (k - 1) • (k - 2)/2 + M (13) 
q(j, k, M) = j + (k - 1) • (k - 2)/2 + M(M + l)/2 (14) 
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by choosing j and k according to the condition: 

j<k k = 2...M (15) 

Differently from the case of the planar beam combiner, the matrix for the DBC has in 
general non-zero coefficients. A sufficient condition to use of the array as an interferometric 
beam combiner is that the corresponding a-matrix is invertible and well conditioned, i.e. 
operations involving the multiplications of the matrix or its inverse have a minimal impact 
on the accuracy of the result. A practical way of gauging the performance of the combiner is 
to evaluate the condition number. For an invertible matrix, the condition number is defined 
as the ratio of the maximum to the minimal eigenvalues. In this sense, the best configuration 
for the combiner is the one featuring the minimal condition number. For a given array, a few 
input configuration can fulfill this requirement for a limited range of propagation distances. 
For combiners featuring a few channels only (< 10), the best configuration can be found by 
means of direct numerical search. 



C. Comparison with ideal combiners 

Numerical simulations of hypothetical interferometric observations (such as those pre- 
sented in [l2, Q) can be used to asses the advantage of one scheme over the other. The 
simulated observations included a photon shot-noise error source and were found to deliver 
performance comparable to existing beam combiners. Here a different approach was followed 
and a comparison of the condition number of the DBC with that of ABCD beam combiners 
was performed. Even though the condition number describes the noise amplification in the 
worst possible scenario only, it has the advantage of being easy to compute and to represent 
their performance with a single number. In Table 1, the DBC with N ranging from 3 to 6 
is compared to an equivalent ABCD combiner. 

The DBC combiner model uses a square lattice geometry and takes into account the 
nearest-neighbor site coupling with unitary strength in the horizontal and vertical direction, 
as well as a diagonal coupling to the next-nearest-neighbor sites with relative strength equal 
to 0.13. This model has been found to describe accurately square arrays of waveguides 
manufactured with the direct laser- writing technique [6]. The optimal input configuration 
and length of the DBCs were found by direct numerical test of all possible configurations. 
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N Baselines kabcd ^dbc Configuration 

3 3 2.0 8.3 

4 6 1.7 7.7 

5 10 2.6 15.0 

6 15 2.2 20.3 

TABLE I. Theoretical performance of the DBC as compared to the ideal ABCD combiners. The 
input configuration for the best DBC is also shown (black waveguides: excitation points). 

Table 1 shows that the condition number of DBCs grows nearly linearly with the number of 
combined telescopes, with the exception of the case N=4 which features a better stability 
than the N=3 combiner. This fact suggests that matching the symmetry of the combination 
problem with that of the array could result in beam combiners of improved performance. 
Indeed, preliminary investigation, has shown that simple deformations of the array (e.g. 
different site separation in the vertical and horizontal direction) can sensibly reduce the 
condition number of the best input configuration of the DBC. Moreover, for the 6-fold- 
DBC, the condition number of the best configuration is smaller for hexagonal lattices. Full 
account of a more systematic investigation of the geometric effects will be given elsewhere. 

In all cases, the ABCD approach provides better conditioning than DBCs. This may 
be related to the fact that the number of interferometric measurements is larger in ABCD 
combiners than in DBC settings, thus resulting in higher stability. To prove this trend, a 
simple test was made by estimating the condition number for over-sized DBC arrays, i.e. 
arrays where the number of waveguides exceeds A^ 2 . For this for the case of ABCD 

combiners, the pseudo-inverse of the a-matrix can be calculated. Table 2 resumes the result 
for the case of a combiner of A^-beams on a square array of (N + l) 2 waveguides. The cases 
N = 3 and N = 4 are examined and show a condition number significantly lower than the 
classical DBC design and slightly higher than the one expected for the ABCD combiner. 

IV. CONCLUSIONS AND PERSPECTIVES 

A few conclusions can be drawn from the result of comparison outlined in the previous 
section. A first conclusion is that ABCD combiners are better conditioned than DBCs. The 
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N Array kdbc Configuration 

•odd 

•ooo 

3 4x4 3.0 8888 

••000 

OOOO* 

OOOOO 

4 5x5 4.3 888S8 

TABLE II. Examples of oversized DBC arrays. The condition number in this case is closer to the 
values expected for the ABCD combiners (See Table 1). 

price for the better performance is that the design of the component become increasingly 
difficult as number of baselines and combined telescopes grows. Concerning DBCs, original 
arrangements for N— telescopes combination featuring A^ 2 arrays have condition number 
growing linearly with the number of combined beams, which makes them uncompetitive to 
ABCD (in terms of the condition number). However, DBCs featuring M > N 2 waveguides 
show a dramatic reduction of the condition number making them comparable to the ABCD 
combiners. 

Considering these results, it is possible to identify the potential of the DBC approach 
along three main lines, i.e. sensitivity, fabrication and flexibility. Assuming that coupling 
losses in the structure are negligible, the DBC shares with the planar ABCD combiners a 
potentially higher sensitivity as compared to methods based on spatial multiaxial combiners 
such as in AMBER [16] or in pupil remapping instruments [lOj, because the number of pixels 
required to read out the coherence information is much smaller. The ease of fabrication of 
DBC arrays bears the greatest potential advantage of the DBC over other methods. In 
fact, as opposed to intricate networks of couplers and cross-overs, DBCs are simple arrays of 



waveguides which could be fabricated by rod-in-tube multicore fiber draw 19| or direct laser 

on 

writing p, 122|. The last technique entails also the potential flexibility of the DBC approach, 
since direct laser writing can be used to write photonic components on a large variety of 



materials including fused silica 15| and chalcogenide glasses [23|, thus offering the possibility 
to cover a wide spectral band from visible to mid infrared. Especially interesting would be 
the realization of such components for interferometric imaging in the mid infrared, with 
obvious applications to the astrophysics of sub-stellar companions or protoplanetary disks. 
Disadvantages of this technique are that the uniformity of the waveguides may critically 
depend on the stability of the femtosecond laser source during the writing process. Initial 
experiments however have shown that a non-uniformity of the coupling strength across the 
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array of ±10% do not irremediably compromise the performance of a real DBC 14(. Work 
is currently in progress to improve the results of laboratory tests, in view of a real on sky 
demonstration of the capabilities of the DBC concept. 
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